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ABSTRACT 

Wc use the corrections An? to the standard Newtonian/Einstcinian pcrihchon preces- 
sions of the inner planets of the solar system, recently estimated by E.V. Pitjeva by 
fitting a huge planetary data set with the dynamical models of the EPM ephemerides, 
to put constraints on the position of a putative, yet undiscovered large body X of mass 
Mx, not modelled in the EPM software. The direct action of X on the inner planets 
can be approximated by a Hooke-type radial acceleration plus a term of compara- 
ble magnitude having a fixed direction in space pointing towards X. The perihelion 
precessions induced by them can be analytically worked out only for some particu- 
lar positions of X in the sky; in general, numerical calculations are used. We show 
that the indirect effects of X on the inner planets through its action on the outer 
ones can be neglected, given the present-day level of accuracy in knowing A-dj. As 
a result, we find that Mars yields the tightest constraints, with the tidal parameter 
/Cx = GMx/rx < 3 X 10"^'' s~^. To constrain rx we consider the case of a rock-ice 
planet with the mass of Mars and the Earth, a giant planet with the mass of Jupiter, 
a brown dwarf with Afx = SOmjupitor, a red dwarf with M = O.SM© and a Sun- 
mass body. For each of them we plot rx'" as a function of the heliocentric latitude (3 
and longitude A. We also determine the forbidden spatial region for X by plotting its 
boundary surface in the three-dimensional space: it shows significant departures from 
spherical symmetry. A Mars-sized body can be found at no less than 70 — 85 AU: such 
bounds are 147 - 175 AU, 1, 006 - 1, 200 AU, 4, 334 - 5, 170 AU, 8, 113 - 9, 524 AU, 
10, 222 — 12, 000 AU for a body with a mass equal to that of the Earth, Jupiter, a 
brown dwarf, red dwarf and the Sun, respectively. 

Key words: Solar system objects; Low luminosity stars, subdwarfs, and brown 
dwarfs; Kuiper belt, trans-Neptunian objects; Oort cloud; Celestial mechanics 



1 INTRODUCTION 

Does the solar system contain an undiscovered massive 
planet or a distant stellar companion of the Sun? 

The history of an hypothetical Planet X dates back to 
the early suggestions by the astronomer P. iLowelll (|l915l ) 
who thought that some glitches in the orbit of Uranus 
might be caused by what he dubbed Planet X. In 1930, 
the search that L owell initiated led to the discovery of Pluto 
l|Tombaughlll96ll ). The various constraints on the mass and 
position of a putative Planet X, as candidat e to accommo - 
date the alleged orbit al anomalies of Ur anus (|Bruninill 19921 ). 
were summarized bv lHogg et al.l l|l99ll ). Later, th e claimed 
residu als in the orbit of Uranus were explained bv lStandishI 
l|l993t ) in terms of small systematic errors as an underes- 
timate of the mass of Neptune by 0.5%. However, the ap- 
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peal of a yet undiscovered solar system 's body of planetary 
size n ever faded. Indeed, according to iLvkawka fc Tadashil 
(|2008l ). a ninth planet as large as the Earth may exist be- 
yond Pluto, at about 100-170 AU, to expla in the architecture 
of the Edgeworth-Kuiper Belt. Previouslv. iBrunini fc Melital 
(2002) proposed the existence of a Mars-size body at 60 
AU to explain certain features in the distribution of the 
Trans-Neptunian Objects (TNOs) like the so-called Kuiper 
Cliff where low-eccentricity and low-inclination) Kuiper Belt 
objects (KBOs) with semimajor axes greater than 50 AU 
rapidly falls to zero, although several problems in explain- 
ing other features of the Kuiper Belt with su ch a hypothe- 
sis were pointed out later l|Melita et al.|[20o3 ). According to 
iMatese et all l| 19991 ). a perturber body of mass m « 1.5mjup 
at 25 kAU would be able to explain the anomalous distri- 
bution of orbital elements of approximately 25% of the 82 
new class I Oort cloud com ets. A similar hypothesis was 
put forth by iMurravl (|l999h : for a more skeptical 



view, see 
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iHorner fc EvansI (|2002l ). iGomes et all (120061 ) suggested that 
distant detached objects among the TNOs (perihehon dis- 
tance g > 40 AU and semimajor axis, a > 50 AU) may have 
been generated by a hypothetical Neptune-mass companion 
having semiminor axis b < 2, 000 AU or a Jupiter-mass com- 
panion with & < 5, 000 AU on significantly inclined orbits. 
However, generally speaking, we stress that we still have 
bad statistics about t he TNOs and the E dgeworth-Kuiper 
Belt objects (.Sebastien fc Morbidellil l2007l : iSchwamb et al.i 
I2OO9I ). 



Concerning the existence of a putative stellar compan- 
ion of the Sun, it was argued with different approaches. As 
an explanation of the peculiar properties of certain pulsars 
with anomalously small period derivatives, it was suggested 
by [Harrison (1977) that the barycenter of the solar system 
is accelerated, possibly because the Sun is a member of a 
binary system and has a hitherto undetected comp anion 
star. Latest studies by IZakamska fc Tremaind (|2005l ) con- 
strain such a putative acceleration at a Abary ^ 1-2 x 10~^ 
m s~^ level (4 x 10~^* s~^ in units of ^barv/c, where c is 
the speed o f light in vacuum) . IWhitmire fc JacksonI (Il984l ) 
and I David (|l984f ) suggested that the statistical periodic- 
ity of about 26 Myr in extin ction rates on the Earth over 
the last 250 Myr reported by iRaup fc Scpkoskil (|l984f ) can 
be explained by a yet undetected companion star (called 
Nemesis) of the Sun in a highly elliptical orbit that peri- 
odically would disturb comets in the Oort cloud, causing 
a large increase in the number of comets visiting the in- 
ner solar system with a conseque ntial increase in impact 
events on Earth. In a recent work, iMullerl (|2002h used the 
hypothesis of Nemesis to explain the measurements of the 
ages of 155 lunar spherules from the Apollo 14 site. The 
exact nature of Ne mesis is uncertain; it could be a red 
dwarf (iMullej |2002|) f0.075 ^ m < 0.5 Mq), or a brown 
dwarf (jWhitmire fc JacksonI Il984l ) (m « 75 - SOmjup). 
For some resonant mechanisms between Nemesis and the 
Sun triggering Oo rt cloud's com et showers at every peri- 
helion passage, se e Vandervoprt fc Sather (19 9311. Recently, 
iFoot fc Silagadzd (120011 ) and lSilagadzd (|200ll ) put forth the 
hypothesis that Nemesis could be made up of the so called 
mirror matter, whose existence is predicted if parity is an 
unbroken symmetry of nature. 



In this paper we constrain the distance of a very distant 
body for different values of its mass in a dynamical, model- 
independent way by looking at the gravitational effects in- 
duced by it on the motions of the inner planets orbiting in 
th e 0.4 — 1.5 AU range. The s ame approach was follow ed 
bv lKhriplovich fc Pities! (120061 ) and iKhriplovichI (|2007l ) to 
put constraints on density of diffuse dark matter in the solar 
system. In Section [5] we calculate the acceleration imparted 
by a distant body X on an inner planet P and the resulting 
perihelion precession averaged over one orbital revolution 
of P. In Section |3] we discuss the position-dependent con- 
straints on the minimum distance at which X can exist for 
several values of its mass, and depict the forbidden regions 
for it in the three-dimensional space. In Section [3] we com- 
pare our results to other constraints existing in literature 
and summarize our findings. 



2 THE PERIHELION PRECESSIONS 

INDUCED BY A DISTANT MASSIVE BODY 

The gravitational acceleration imparted by a body of mass 
Mx on a planet P is, with respect to some inertial frame, 

GMx 



^px = 



rx - rp 



-3 (rx-rp). 



(1) 



If, as in our case, it is supposed rx 2> rp, by neglecting terms 
of order 0(rp/rx), it is possible to use the approximated 
formula 



1 



rx - rp 



1 + 



3 (rp • nx) 



where 



nx 



rx 



rx 

rx 



+ 



(2) 



(3) 



is the unit vector of X which can be assumed constant over 
one orbital revolution of P. By inserting eq. (|2)| in eq. ([1} 
and neglecting the resulting term of order 0{rp/r^) one has 

. GMx 3GAfx (rp • nx) - GMx. , 
Apx ~ 5— H q nx H 5 — nx. (4) 

rl rl rl 

The acceleration of eq. (|4]) consists of three terms: a non- 
constant radial Hooke-type term, a non-constant term di- 
rected along the fixed direction of fix and a constant (over 
the typical timescale of P) term having the direction of fix 
as well. An equation identical to eq. @ can be written also 
for the Sun by replacing everywhere rp with rs; thus, since 
we are interested in the motion of the planet P relative 
to the Sun, the constant term cancels out, and by posing 
r = rp — rs one can writes down the heliocentric perturb- 
ing acceleration felt by the planet P due to X 

GMx 3GMx (r • fix) . 
Ax = Anooko + ^nx = 3 — rH 3 nx. (5) 

^x ''x 

In fact, the planetary observation-based quantities we will 
use in the following have been determined in the frame of 
reference of the presently known solar system's baricenter 
(SSB). By the way, eq. ((SJ, where r is to be intended as the 
position vector of P with respect to the known SSB, is valid 
also in this case. Indeed, by repeating the same reasonings 
as before, we are led to a three-terms equation like eq. ^ 
in which the vectors entering it refer to the SSB frame; now, 
if X existed, the SSB frame would be uniformly accelerated 
by 

A GMx . 

AssB = — 2 — "X, (6) 
''x 

so that P, referred to such a non-inertial SSB frame, would 
also be acted upon by an inertial acceleration 



GMx . 
9— "■x, 



(7) 



in addition to the gravitational one of eq. ([4]), which just 
cancels out the third constant term in eq. ^ leaving us 
with eq. Note that, in general, it is not legitimate to 
neglect A„^ with respect to AhooIcc since 

An 



— 3cosi/), 



(8) 



^Hookc 

where tp is the angle between r and nx which, in general, 
may vary within < tp < 2n during an orbital revolution 
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of P. Moreover, if the acceleration of a distant body could 
only be expressed by a Hooke-type term, this would lead to 
the absurd conclusion that no bodies at all exist outside the 
solar system because the existence of an anomalous Hooke- 
like acceleration has been ruled out b y taking th e ratio of 
the perihelia of different pairs planets l|loridl2008l l. 

The acceleration of eq. ([5]) is in agreement with the 
potential 



Ux = 



GMx 
2ri 



3 (r • hx) 



(9) 



proposed bv lHogg et all l|l99ll ). 

The action of eq. ^ on the orbit of a known planet 
of the solar systet ri can he treated pe rturbatively with the 
Gauss equations jBertotti et alj 120031 ) of the variations of 
the Keplerian orbital elements 



(10) 



sin / + At ^aasf + " (^1 " 



cos / + 1 + - sin / - cos J 



. —aJ-] - -Jl - c2 — +COSI — 
Tia V a / \ dt dt 



where a, e, I, Q, ui and A4 are the semi-major axis, the 
eccentricity, the inclination, the longitude of the ascending 
node, the argument of pericenter and the mean anomaly of 
the orbit of the test particle, respectively, u — lj + f is the 
argument of latitude (/ is the true anomaly), p — a(l — 
e^) is the semi-latus rectum and n = ^ GM/a^ is the un- 
perturbed Keplerian mean motion; Ar,At, A„ are the radial, 
transverse and normal components of a generic perturbing 
acceleration A. By evaluating them onto the un-perturbed 
Keplerian ellipse 

a{l-e^) 



■ l + ecos/ ^^^^ 

and averaging] them over one orbital period Pb of the planet 
by means of 



dt 



(1 



27r(l -f ecos/)- 



:df, 



(17) 



it is possible to obtain the long-period effects induced by eq. 

In order to make contact with the latest observational 
determinations, we are interested in the averaged rate of 
the longitude of the pericenter = cij + O; its variational 
equation is 



^f + At 



n/ + 2s 



Indeed, the astronomer iPitieval (|2p05') has recently es- 
timated, in a least-square sense, corrections Azu to the 

^ We used u = f because over one orbital revolution the pericen- 
tre u! can be assumed constant. 



standard Newtonian/Einsteinian averaged precessions of the 
perihelia of the inner planets of the solar system, shown in 
Table [T] by fitting almost one century of planetary observa- 
tions of several kinds with the dynamical force models of the 
EPM ephemerides; since they do not include the force im- 
parted by a distant companion of the Sun, such corrections 
are, in principle, suitable to constrain the unmodelled action 
of such a putative body accounting for its direct action on 
the inner planets themselves and, in principle, the indirect 
one on them through the outer planetqj; we will discuss such 
an issue later. 

Concerning the calculation of the perihelion precession 
induced by eq. (O, as already noted, rx can certainly be 
considered as constant over the orbital periods Pb ^ 1.5 yr 
of the inner planets. Such an approximation simplifies our 
calculations. In the case of the Hook-type term of eq. ((S} 



-JCxr, tCx 



GMx 



(19) 



the task of working out its secular orbital effects has been 
already performed several times in literature; see, e.g.. Iloriol 
(.2008. ) and references therein. The result for the longitude 
of perihelion is 



<^)hoo 



3r- y]^^ 
2 n 



3 GMx 
2 rl 



GMa 



(20) 

The treatment of the second term A^^ of eq. (O is much 
more complex because it is not simply radial. Indeed, its 
direction is given by fix which is fixed in the inertial space 
{x,y,z} during the orbital motion of P. Its components are 
the three direction cosines, so that the cartesian components 
of 



Axi -\- Ayj ~\- Azh 



are 
Ax 



2>K.x{xnx -\- yuy + znz)nx 



SlCxixux + yUy + znz)ny 



(21) 



(22) 



(23) 



Az = SICxixnx + yUy + zn2)n2. (24) 

In terms of the standard heliocentric angular coordinate^ 
A, l3, we can write the components of fix as 



1 — sin^ Px cos Ax , 



\l 1 ~ sin^ /3x sin Ax , 



Uz = sin f3x , 



(25) 
(26) 
(27) 



with nx+Uy+Tiz = 1. In order to have the radial, transverse 



^ Only their standard Af— body mutual interactions have been, 
indeed, modelled in the EPM ephemerides. 

^ Recall that A is the usual longitude of the standard spherical 
coordinates, while f3 = 90° — 9, where 9 is the usual co-latitude 
of the standard spherical coordinates. 
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and normal components of A„y_ to be inserted into the right- 
hand-sides of the Gauss equations 



At = A„ 



Afi — An 



(28) 



(29) 



(30) 



we need the cartesian components of the co-movin g unit 
vectors r,i,n: they are (jMontenbruck and Gilll[2000t ) 



cos f2 cos u — cos / sin f2 sin u 
sin Q cos u + cos I cos Q sin u 
sin I sin u 

— sin It cos fl — cos / sin Q cos it 

— sin SI sin u -I- cos I cos f2 cos it 

sin 7 cos It 

(sin / sin Q \ 
— sin / cos Q I 
cos / / 



(31) 



(32) 



(33) 



In our specific case, Q is the longitude of the ascending 
node which yields the position of the line of the nodes, i.e. 
the intersection of the orbital plane with the mean ecliptic at 
the epoch (J2000), with respect to the reference x axis point- 
ing toward the Aries point T. Since, according to eq. (|22p -ea. 
H24p . the planet's coordinates x,y,z enter the components 
of A„y^, we also need the expressions for the unpertu rbed 
coordinates of the planet (jMontenbruck and Gillll2000l ) 



X — r (cos cos u — cos I sin fl sin u) . 



y — r (sin Q cos u + cos / cos fl sin it) , 



(34) 



(35) 



for its imperfect modelling in the EPM ephemerides. Indeed, 
it is easy to show that the secular precession of the longitude 
of pericenter of a test body moving along an orbit with small 
eccentricity around an oblate body of equatorial radius R is 



J2 



R 



J2 



(1 - e2 



5 2 r 
-cos /- 



cos / 



(38) 

Now, let us consider the case in which X is directed along 
the z axis and the planet P has 1 = 0; only the retrograde 
Hooke-type precession of eq. (|20|l would be non-vanishing, 
while eq. (|38|) would induce a prograde rate. 



3 THE THREE-DIMENSIONAL SPATIAL 

CONSTRAINTS ON THE LOCATION OF X 

In Figure [1] we plot the maximum value of ICx, obtained 
from the perihelion of Mars which turns out to yield the 
most effective constraints, as a function of /3 and A. It turns 
out that /Cx'"' is of the order of 10"^"* s'^ 

The minimum distance rx can be obtained for different 
values of Mx as a function of l3 and A as well according to 



mm 



f GMx \ 



1/3 



(39) 



We will consider the case of two rock-ice bodies with the 
masses of Mars and Earth, a giant planet with the mass of 
Jupiter, a brown dwarf with Mx ~ SOmjupiter, a red dwarf 
with Mx = 0.5mQ and a Sun-mass body with Mx = Mq, 
non necessarily an active main-sequence star. 

It is also possible to visualize the forbidden region for 
X in the three-dimensional space by plotting its delimiting 
surface whose parameteric equations are 



rx'"(/3, A) cos 13 cos A, 



(40) 



2 = r sin J sin u. 



(36) 



Thus, Ar,At,A„ are non- linear functions of the three un- 
known parameters /Cx,Ax,/3x of X; by computing the av- 
eraged perihelion precessions induced by ^Hooke and An^ 
and comparing them with the estimated corrections Ari; to 
the usual perihelion precessions it is possible to have an up- 
per bound on /Cx, and, thus, a lower bound on rx for given 
values of Mx, as a function of the position of X in the sky, 
i.e. rx'" = fx'"(Ax, /3x). Relatively simple analytical expres- 
sions can be found only for particular positions of the body 
X; for example, by assuming rix = Uy = 0, i.e. by consider- 
ing X located somewhere along the z axis, it is possible to 
obtain 



g/CxVl-e^sin^J 
4n 



^ cos 2lo 



(37) 



In Table [5] we quote the lower bound on rx = zx for several 
values of Mx obtained from the maximum value of the es- 
timated correction to the standard rate of the perihelion of 
Mars applied to eq. (|20p and eq. (|37p . In general, the prob- 
lem must be tackled numerically; also in this case it turns 
out that the data from Mars yield the tightest constraints. 
However, eq. (|20|l and eq. (|37p are already enough to show 
that the effect of a body X could not be mimicked by a cor- 
rection to the Sun's quadrupole mass moment J2 accounting 



rx'"(/3. A) cos (3 sin A, 



rx'"(/3,A) sin/3. 



(41) 



(42) 



In Figure [5] we plot the minimum distance at which a 
body of mass Mx = rnMa.rs can be found as a function of its 
latitude and longitude (3 and A. The largest value is about 
85 AU and occurs in the ecliptic plane {f3 = 0) at about 
A = 0, 60, 150, 250, 345 deg. Note that for (3 = ±90 deg, i.e. 
for X along the z axis, the minimum distance is as in Ta- 
ble [5J i.e. 70 AU. For just a few positions in the sky such a 
Mars-sized body could be at no less than 20 AU. In Figure 
Owe depict the surface of minimum distance delimiting the 
spatial region in which such a body can exist according to 
the data from Mars. It has not a simple spherical shape, as 
it would have had if X exerted an isotropic force on Mars, 
and it has a precise spatial orientation with respect to the 
{x,y,z} frame. The analytical results of Table [2] obtained 
with eq. (|37p in the case /3 = ±90 deg are confirmed by 
Figure [5] and Figure O By slicing the surface with a vertical 
symmetry plane, it turns out that the shape of the central 
bun is rather oblate, approximately by « 0.77, contrary to 
the lateral lobes which are more spherical. Moreover, the 
largest forbidden regions in the ecliptic plane are two ap- 
proximately orthogonal strips with a length of about 180 AU 
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(see Figure |4] depicting the situation in the ecliptic plane), 
while in the vertical direction there is a strip approximately 
120 AU long. 

Figures analogous to Figure [2]-Figure |4] hold for bodies 
with the mass of the Earth, Jupiter, a brown dwarf (m = 
SOmjupiter), a red dwarf with M = 0.5Mq and a Sun-sized 
body. The results of Table [2] concerning the position of X 
along the z axis are confirmed. The same qualitative features 
of the case Mx = mMars occur. Concerning an Earth-sized 
body, it could mainly be found at no less than 147—175 AU, 
with a minimum distance of 35 AU for just a few positions 
in the sky. The minimum distance of a Jupiter-like mass is 
1, 006 - 1, 200 AU, with about 200 AU in some points. For 
a brown dwarf (Afx ~ 80 mjupitor) the limiting distance 
is mainly 4, 334 — 5, 170 AU, with a minimum of 861 AU 
at some positions, while for a red dwarf {Mx ~ 0.5 A/0) 
it is 8, 113 - 9,524 AU, with a lowest value of 2,000 AU. 
Finally, a Sun-mass body cannot be located at less than 
10, 222-12, 000 AU for most of the sky positions, with 2, 520 
AU at just a few points. In the case of an Earth-sized body, 
the length of the largest ecliptical forbidden strip turns out 
to be about 400 AU, while the vertical one amounts to about 
300 AU. 

Concerning the strategy adopted so far, let us note that 
we compared the estimated Aw for the inner planets to their 
computed perihelion precessions due to X through eq. ([5]); 
in fact, one should, in principle, also take into account the 
indirect effects of X on the inner planets of the solar system 
through its direct action on the outer planets. Indeed, if one 
of them is perturbed by X, its action on the inner planets will 
differ from the standard A'^— body one, fully modelled in the 
EPM ephemerides. To roughly evaluate such an efi'ect, let 
us reason as follows. The position of both X and of an outer 
planet acted upon by X can approximately be considered 
as constant in space with respect to a rocky planet over a 
typical orbital period Pb ^ 1 yr. The maximum value of the 
disturbing acceleration imparted by X on a giant planet G 
like Jupiter will be of the order of 



7 X 10"^^ m s"^ 



(43) 



we used JCx = 3 x 10 s ^. Aqx will displace the position 



of G with respect to its usual position by an amount 

Agx 



Arc 



3 X 10'' m 



(44) 



over f'b 1 yr during which Aqx can be considered as 
uniform. As a consequence, the maximum extra-acceleration 
Ay4pG induced on an inner planet P by G during Pb can be 
evaluated as 



AA 



PG 



12 



Arc 1 r < 2 X lO"^'' m s"^ 



(45) 



we have just used eq. ((5]) adapted to the present case. The 
perihelion precession induced by AApc can roughly be eval- 
uated by dividing it by the orbital velocity of P, i.e. an, so 
that 

, . , AApG 



< 10 ^ arcsec cty ^. 



(46) 



It is about one order of magnitude smaller than the present- 
day accuracy in knowing the extra-precessions of the peri- 
helia of the inner planets, as shown by Table [1] Thus, we 
conclude that the indirect effects of X on the inner planets 



through its action on the giant ones can be neglected, as we 
did. 



4 DISCUSSION AND CONCLUSIONS 

We will, now, compare our dynamical constraints with other 
ones obtained with different methods. 

First, we compare the dynamical constraints of Table [2] 
and, more generally, of our full analysis with those obtain- 
able as rx = \/ GA'fx /^bary from the upper bound on the 
solar system barycenter' s acceleration Abarv < 1 . 2 x 10 ~^ m 
s~^ recently derived bv IZakamska &: Tremain"3 (|2005l ) with 
an analysis of the timing data of several millisecond pulsars, 
pulsars in binary systems and pulsating white dwarf. Indeed, 
an acceleration exerted on the known SSB would affect the 
observed value of the rate of the period change of astronom- 
ical clocks such as pulsars and pulsating white dwarfs. Such 
a method is able to provide a uniform sensitivity over the 
entire sky. The pulsar-based, isotropic 1/rx constraints are 
summarized in Table [3] With "isotropic" we mean that the 
value of A^^ used can be ruled out for 100% of the sky at 
95 % confidence level with practically all the methods used 
by Zakamska fc Tremai"n3 (120051 '): by using PSR B1913 -I- 6 
the quoted acceleration is ruled out at 95% confidence for 
94% of the sky. They are not competitive with those of Table 
[2] and of the rest of our an alysis fo r all t he bodies consid- 
ered. However, according to lEichleil (|2009l ). future timing of 
millisecond pulsars looking for higher order pulsar period 
derivatives should be able to extend such limits to several 
thousand AU within a decade. A precise limit on the unmod- 
eled relative acceleration between the solar system and PSR 
J0437-4715 has recently been obtained byJDeller (2008) by 
comparing a VLBI-based measurement of the trigonometric 
parallax of PSR J0437-4715 to the kinematic distance ob- 
tained from pulsar timing, which is calculated from the pul- 
sar's proper motion and apparent rate of change of orbital 
period. As a result, Jupiter-mass planets within 226 AU of 
the Sun in 50% of the sky (95% confidence) are excluded. 
Proposals to search for primordial black holes (PBHs) with 
the Square Kilometer Array by using modification of pulsar 
timing residuals when PBHs pass within about 1,000 AU 
and impart impulse acceler a tions to the Earth have been 
put forth bv lSeto fc Cooravl |20o3). 

Other acceleration- type 1/rx constraints were obtained 
by looking for a direct action of X on the Pioneer 
10/1 1 spacecraft from an inspection of their tracking data 
(| Anderson , 198&1 : the upper bound in the SSB acceleration 
obtained in this way is 4.2 x 10"^" m s~^. This translate in 
a re-scaling of the values of Table [3] by a factor 1.7, still not 
com petitive wit h ours. 

iHogg et al.] (119911 ) looked for tidal- type l/r'i con- 
straints, like ours, with a dynamical approach based on nu- 
merical simulations of the data of the four outer planets over 
the time span 1910—1990. They found a relation among Afx, 
in units of terrestrial masses, rx, in units of 0.1 kAU, and 
(70, which is the standard deviation of the assumed Gaussian 
random errors in A and /3, in units of 0.1 arcsec; it is, for the 
ecliptic plane. 



( ) 
VlOO AU/ 



A/x \ / 0.1 arcsec 



-,1/3 



(6 - 10) 



(47) 



6 L. lorio 



By assuming, rather optimistically, cro ~ 0.01 arcsec, eq. 
(|47)l yields for Mx = m® rx > 120 AU, which is generally 
smaller than our limits for the ecliptic plane; a more conser- 
vative value of (To ~ 0.1 arcsec yie l ds rx > 55 AU. It may be 
interesting to recall that iBruninil (|l992l ) predicted the exis- 
tence of an alleged perturber of Uranus with Mx = 0.6m^ 
at a heliocentric distance of 44 AU. 

Let us, now, focus on direct observational limits. A 
planetary body would reflect the solar light and, therefore, 
could be detecte d in the optical or near-infrared surveys. 
iTombaughl l|l96ll ) conducted an all-sky optical survey that, 
in the plane of the ecliptic, yielded the following constrain 

(lOO Au) ~ .(^) 3A ' ^^^^ 

by assuming a visual albedo p = 0.02 (|Hogg et al ] ll99ir) . For 
an Earth-sized body the limit is rx > 81 AU. Anoth er opti- 
cal su rvey yielding similar results was performed bv iKowaJ 
In addition to the reflected sunlight, also the X's 
own thermal emission would be detecta ble in the mid-to 
far-in frared. The IRAS All-Sky Survey (jNeugebauer et al] 
would have been able to discover a gas giant in the 
range 70 AU< rx < 400 AU, but without success, in agree- 
ment with our bounds for Jupiter -sized bodies. A recent 
ecliptic survey was done bv iLarsenl (j2007) with the Space- 
watch telescope. This survey was sensitive to Mars-sized ob- 
jects out to 300 AU and Jupiter-sized planets out to 1,200 
AU; for low inclinations to the ecliptic, it ruled out more 
than one to two Pluto-sized objects out to 100 AU and one 
to two Mars-sized objects to 200 AU. Concerning the case 
of a Sun-mass body, the presence of a main-sequence star 
above the hydrogen-burning limit was e xcluded within 1 pc 
by the all-sky synoptic Tycho-2 survey (|H0g et al.ll2000l ). 

Observational constraints on the properties of X might 
be obtained by sampling the astrometric position of back- 
ground stars o ver the entire sky with the future astrometric 
GAIA mission (jScott et al.ll2005h . Indeed, the apparent mo- 
tion of X along the parallactic ellipse would deflect the angu- 
lar position of distant stars due to the astrometric microlens- 
ing ("induced parallax"). A Jupiter-sized planet at 2,000 AU 
in the ecliptic plane could be detected by G AIA. Also the 
phenomenon of mesolensing (|Di Stefandl200^ could be used 
to gain informatio n on p ossible planets at distances > 1, 000 
AU. lBabich et al I ("2007) proposed to use the spectral distor- 
tion induced on the Cosmic Microwave Background (CMB) 
by putative distant masses to put constraints on the physical 
properties and distanc es of them. W ith the all-sky synoptic 
survey Pan-STARRS (|jewittll2003l ) massive planets such as 
Neptune would be undetectable through reflected sunlight 
beyond about 800 AU, while a body with Mx = O.IM© 
would be undetectable for rx > 2, 000 AU. 

In conclusion, the dynamical constraints on a still undis- 
covered planet X in the outer regions of the solar system we 
dynamically obtained from the orbital motions of the inner 
planets of the solar system with the extra-precessions of the 
perihelia estimated by Pitjeva with the EPM ephemerides 
are tighter than those obtained from pulsar timing data 
analysis and the outer planets dynamics; moreover, they 
could be used in conjunction with the future planned all- 
sky surveys in the choice of the areas of the sky to be in- 
vestigated. If and when other teams of astronomers will in- 
dependently estimate their own correction to the standard 



perihelion precessions with different ephemerides it will be 
possible to repeat such tests. Moreover, a complementary 
approach to that presented here which could also be imple- 
mented consists of modifying the dynamical force models 
of the ephemerides by also including the action of X on all 
the planets of the solar system, and repeating the global flt- 
ting procedure of the entire planetary data set estimating, 
among other parameters, also those which directly account 
for X itself and looking at a new set of planetary residuals. 
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Table 1. Estimated corrections Aro, in milliarcsec cty~^ (1 arc- 
sec cty~^ = 1.5 X 10~^^ to the standard Newton/Einstein 
pe rihelion precessions of the inner planets according to Table 3 
of iPitieval ||2005|') (Mercury, Earth, Mars). The result for Venus 
has been obtained by recently processing radiometric data from 
Magellan spacecraft (E.V. Pitjeva, private communication, 2008). 
The errors are not the formal, statistical ones. The SSB frame, as- 
sumed as inertial, i.e. without modeUing the action of a putative 
body X, has been used. 



Mercury 


Venus 


Earth 


Mars 


-3.6 ±5.0 


-0.4 ±0.5 


-0.2 ±0.4 


0.1 ±0.5 



Table 2. Minimum heliocentric distance rx, in AU, at which a 
still unseen object having a mass Mx equal to that of the as- 
tronomical bodies listed here can be located along the z axis 
(/3 = ±90 deg) according to eq. I|20| l. eq. I|37| l and the maximum 
value of the extra-precession of the perihehon of Mars, according 
to Table [T] 



Mars 


Earth 


Jupiter 


Mx = 80mjup 


Mx = ^ 


Sun 


70 


147 


1,006 


4,336 


8, 113 


10, 222 



Table 3. Approximate minimum hehocentric distance rx, in AU, 
at which a still unseen object having a mass Mx equal to that 
of the astronomical bodies listed here can be located accord- 
ing to the hmit Abary = GMx/r^ < 1.2 X 10~^ m s"^ on 
the solar system barycenter' s acceleration recently obtained by 
IZakamska &: Tremainel ||2005| ) from pulsar timing data. We used 
the maximum value A]-,ary/c = 4x 10~^* where c is the speed 
of hght in vacuum, of the SSB acceleration allowed by pulsar tim- 
ing data for about 100% of the sky at 9 5% level of confidence 
(Table 2 bv lZakamska fc Tremain3 ||2005| ')'). 



Mars 


Earth Jupiter 


Sun 


1.4 


3.9 68.6 


2,224 
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Figure 1. Maximum value, in s~^, of ASx = GMx/fx ^ ^ func- 
tion of tlie heliocentric longitude A and latitude of the putative 
body X; its upper bound is of the order of 3 X IQ-^* s-^. The 
perihelion of Mars has been used (Table [TJ . 
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Figure 2. Minimum distance r™'" at which a rock-ice planet 
of mass Mx = niMars can exist as a function of its heliocentric 
longitude A and latitude (3. The perihelion of Mars has been used 
(Table [TJ. 
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X Mars 




y(AU) ^(Ay) 



Figure 3. Surface delimiting the spatial region in which a rock-ice 
planet of mass Mx = mMars can exist according to the dynamical 
constraints from the perihelion of Mars (Table [TJ. The region 
inside the surface is forbidden: the region outside the surface is 
allowed. The red and blue lines correspond to constant values of 
/3 and A, respectively. 
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x(AU) 



Figure 4. Ecliptic view: a Mars-sized rock-ice body can only 
exist outside the region delimited by the red contour. 



